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The universal instability in general geometry 
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The “universal” instability has recently been revived by Landreman, Antonsen and 
Dorlandi, who showed that it indeed exists in plasma geometries with straight (but 
sheared) magnetic field lines. Here it is demonstrated analytically that this instability 
can be present in more general sheared and toroidal geometries. In a torus, the 
universal instability is shown to be closely related to the trapped-electron mode, 
although the trapped-electron drive is usually dominant. However, this drive can be 
weakened or eliminated, as in the case in stellarators with the maximum-J property, 
leaving the parallel Landau resonance to drive a residual mode, which is identified as 
the universal instability. 

PACS numbers: 52.35.Kt, 52.55.Hc, 52.25.Dg 


1 


In an ironic turn of events, the “universal” instability has recently been revivecii. This 
instability of ordinary drift waves was hrst predicted in the 1960but brought into dis¬ 
repute in the late 1970’s by a series of papers demonstrating that it is actually absent in 
the simplest limit of small ion gyroradius, k_i_pi -C 1, if there is no temperature gradient 
or magnetic-held curvature^'-. It was even thought (erroneously) that stability prevails at 
hnite k±pi, so that a plasma without temperature gradients would always be stable in a 
simple sheared magnetic held^. For more information about the long and tortuous history 
of drift-wave (in)stabihty, we refer the reader to the reviews by Horton^ and Connor-. In 
the latest twist of the story, the recent work of Landreman, Antonsen and Borland has now 
shown that the universal instability exists after all, if k±pi is hnite. Their (numerical) proof 
was restricted to a straight, sheared magnetic held, and thus raises the question of whether 
the instability also exists in more general magnetic geometry. 

The universal instability is electrostatic in nature and requires both the ions and the 
electrons to be treated kinetically. The dispersion relation is obtained from the gyrokinetic 
formulation of the quasineutrality condition, which in a hydrogen plasma can be written as 

Sm = SUe, (1) 


where the perturbed density of species a is 


OHa = - - - h / PaJod V. (2) 

The perturbed distribution functions Pa are here dehned such that the full distribution 
functions are fa = /ao(l — Gaf/Ta) + (y'a, fhe electrostatic potential is denoted by 0, the argu¬ 
ment of the Bessel function is k_[V±/VLa (but the electron gyroradius will be neglected), the 
gyrofrequency is denoted by Bq = CaB/ma, and the rest of the notation is standard. We con¬ 
sider the instability in the collisionless, electrostatic approximation, where the distribution 
functions satisfy the gyrokinetic equation 


iV\\Vliga + {UJ- Uda)ga = ^Jq “ (^Ja) faO 
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( 3 ) 


Here, we use the kinetic energy e = mav‘^/2 and A = v‘j_/{v‘^B) as independent velocity- 
space variables, and the parallel derivative Vy is taken holding these constant. The mode 
frequency is denoted by cn, the drift frequency by Uda = k ■ v^a, the diamagnetic frequency 
by u^a = iTa/naGaB‘^)(k± X B) ■ Vria, and we have written 
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with = mav‘^/2Ta and rja = dlnTa/dlnna- The universal instability is most readily ob¬ 
tained in the limit vn/L^^ u, where vxi denotes the ion thermal speed and Ly = l/Vy ln0 
is the length scale of the instability along the magnetic held. We note that parallel ion 
motion can be retained perturbatively, as in previous works, but does not signihcantly alter 
the subsequent analysis. It is thus neglected here, but details of how it affects the argument 
are displayed in the Appendix. For simplicity, we will also neglect Udi', this assumption will 
be discussed more later. The solution to the ion gyrokinetic equation is then 
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and the quasineutrality condition ([T]) becomes 
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where T^ = In{b)e~^ and b = k^pf = k^Ti/{miVt^). 

What remains is the system composed of equation ([3]) for the electrons and the constraint 
(]S|), which is in general difficult to solve analytically. However, the case of a constant, 
unsheared magnetic held is relatively simple to analyze, and we examine this case now as it 
will bring clarity to the general result later. The electron response in this limit becomes 
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if ca = cjj, + ^7 with 7 ~ k\\VTe- It is now a simple matter to hnd unstable regions 


in parameter space from Eq. 
wavelength, kj^pi S> 1, where 


A particularly clear limit is that of small perpendicular 


Tn ~ 


Ti- 
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and thus 
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implying that the plasma is unstable if pi and pe are both less than 2 or both greater than 

2 . 
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We now turn to our main topic of interest, the question of stability in general sheared 
magnetic geometry. For topologically toroidal magnetic helds, the gyrokinetic equation as 
written in Eq. (jS]), can be interpreted in ballooning space. It is worth noting that while 
previous studies of the universal instability have been done in Fourier space, the formula¬ 
tion in ballooning space is equivalent, as can be demonstrated, with some care, by Fourier 
transformation. To pass to ballooning space, we simply write Vy = d/dl, where I denotes 
the arc length along the magnetic held line, now taken to extend inhnitely in both directions. 
As is well known, solutions in the torus are obtained by transformation of the ballooning 
space solutions, but we note that one can also interpret solutions as applying directly to 
inhnite sheared geometries, of which the simplest example is the sheared slab. The following 
analysis is to be interpreted in both ways, i. e. as applying to the old sheared-slab universal 
instability, and also as applying to a generalized toroidal version. 

Because the magnetic held strength and the shear vary along the magnetic held, the 
mode structure and the dispersion relation can no longer be obtained analytically, but we 
can nevertheless extract most of the information we need by the quadratic form obtained 
by multiplying Eq. ([5]) by the complex conjugate of the electrostatic potential, 0*, and 
integrating along the entire held line, —oo<l<oo, 

r [i + F (1 - r„) - (r„ - „6(r„ - r,))] + F r f = o, (7) 

7-00 L ij oj \ B ne 7-00 BJ 

The imaginary part of this equation is 
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where 


/OO di C 

ge(Tv, 

and we conclude that instability is impossible unless Qe{oj)/uj^e is positive. Physically, Qe 
is proportional to the work done by the electrons on the instability^'I^'ii. It was calculated 
in Ref.— and found to be a sum of contributions from passing and trapped particles. 




where, if 0 < 7 Iw^l, the former is given by 
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QepY) =-Tf; - {u - xdx dA ^ |0j(a;,A,a;)| , (9) 
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with 



and 



The corresponding contribntion from the trapped particles is given by 






( 10 ) 


where the snm is taken over all the trapping wells along the held line, denote the locations 
of consecntive bonnce points, and 



We now take the limit oj <C k\\VTei where M -C 1, and note that the passing-particle 
contribntion ([9]), 



( 11 ) 


exactly corresponds to that prodncing the nniversal instability in a straight magnetic held, 
Eq. ([6]), bnt is relatively small compared with the contribntion from trapped electrons. The 
latter is dominated by the m = 0 term in the snm ffTOj) . and becanse of the delta fnnction, 
this term picks ont the resonance u = Hide (an overbar denotes the bonnce average) and is 
larger than Eq. flTT]) by a factor k\]VTel^ 3> 1. 

We are thns led to the conclnsion that in arbitrary toroidal geometry the same driving 
mechanism (parallel Landan resonance) is present to drive an instability as in the case of a 
straight magnetic held, bnt it is generally overwhelmed by the trapped-particle drive. To 
prove that an instability similar to the nniversal one nevertheless exists in certain geometries, 
we consider the limit of steep density gradient, 3> oJde- In this limit, we can neglect ojde 
in a hrst approximation, so that the trapped-electron response simply becomes 



Eqnation d?]) then yields the following variational form for the real freqnencyii, 

uj N[(j)] 


( 12 ) 


uj*e D[(j>y 
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where the functionals N and D are dehned by 


iV[0] = 


,dl 1 rVB. 
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and we have taken = 0 for simplicity. The sums appearing in these expressions are taken 
over all trapping wells along the held line, and 



dl 

v'l - 


denotes the normalized bounce time in such a well. The key idea here is that Eq. (IT^ 
furnishes the zeroth order solution [ur, 0(0)- The growth rate is then computed directly 
from Eq. ([8]), which exactly accounts for the resonances. We do not do this computation 
explicitly, but argue for the existence of a particular solution as follows. 

The denominator D[(j)] in Eq. flT^ is always positive, because of the Schwartz inequality. 
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and so is the numerator N[(j)] when the perpendicular wavelength b = k^pf vanishes. How¬ 
ever, when 6 —)• oo, N[(j)] becomes negative and must therefore vanish at some intermediate 
value of b. There is, therefore, a range of perpendicular wavelengths for which the real 
frequency satisfies 0 < < 1, making the passing-electron contribution to the growth 

rate flTT|) positive if < 2, regardless of the (unknown) mode structure 0(/). As already 
remarked, this contribution is however usually overwhelmed by the trapped-particle drive 
onii. but the latter can be made arbitrarily small by tailoring the magnetic field. In a 
tokamak, this can be achieved by increasing the aspect ratio so as to reduce the fraction of 
trapped particles. In a stellarator, it can be accomplished by minimizing the amount of un¬ 
favorable bounce-averaged curvature. In the limit of a so-called maximum-J stellarator— lii, 
all trapped particles experience average favorable curvature, so that oj^eJXXde < 0 for all or¬ 
bits. Because uj has the same sign as (u^e, there are then no resonances in the delta function 
of Eq. flTOl) and the trapped-electron drive vanishes. In both cases, we are left with an 
instability driven by circulating electrons through a parallel Landau resonance producing 
a response flTT]) very much like that in a straight magnetic held. Note that this instability 
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does not require a temperature gradient but can be driven by a density gradient alone. Fi¬ 
nally, for the purposes of comparing this result with previous work, we note that inclusion 
of finite parallel ion motion in the ion response merely introduces an extra quadratic term 
to the left-hand-side of Eq. ([7]), as shown in the Appendix. This term is of order (k\\VTi/u:Y 
compared with the others, and is additionally small at large kj^pi. 

Let us summarize and discuss our findings. First, we have given new theoretical support 
for the existence of the universal instability in a sheared slab geometry, the context in which 
the mode has been traditionally studied. This may not at once be apparent, but it is because 
the sheared slab is a limit of the gyrokinetic system in ballooning space, and this limit can 
be approached as described above, by continuously deforming a general magnetic field. Our 
argument requires k^pi = 0(1), and so does not contradict stability proofs^^- where this 
parameter was taken to be small. We have also established the existence of the instability 
in topologically toroidal magnetic geometry. In particular, maximum-J stellarators prove 
an intriguing example: It turns out that the stability of trapped-electron modes in these 
conhgurations leaves space for the universal instability to appear. It seems suitably ironic 
that the instability called universal may ultimately be found lurking in the most exotic of 
places. 
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APPENDIX: PARALLEL ION MOTION 


When the first term in Eq. ([3]) is small for the ions, it can be accounted for perturbatively, 
by denoting the lowest-order solution (jl]) by Gi and writing 

g^ = G,- ^V|| (Gi - . 

u \ u J 

Approximating the last term on the right by Gi, we find the ion density perturbation to 
be 

6ni e(j) 6no + 6ni 

rji ^ ^ 
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and in Eq. ([7]) we thus obtain the following additional term on the left-hand side 
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which is a small correction, of order {k\\VTi/ ojY, to the other terms. 
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